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ABSTRACT 

Context. Although debris discs have been detected around a significant number of main-sequence stars, only a few of them are known 
to harbour hot dust in their inner part where terrestrial planets may have formed. Thanks to infrared interferometric observations, it is 
possible to obtain a direct measurement of these regions, which are of prime importance for preparing future exo-Earth characterisa- 
tion missions. 

Aims. We resolve the exozodiacal dust disc around Vega with the help of infrared stellar interferometry and estimate the integrated 
H-band flux originating from the first few AUs of the debris disc. 

Methods. Precise H-band interferometric measurements were obtained on Vega with the 3-telescope IOTA/IONIC interferometer 
(Mount Hopkins, Arizona). Thorough modelling of both interferometric data (squared visibility and closure phase) and spectral en- 
ergy distribution was performed to constrain the nature of the near-infrared excess emission. 

Results. Resolved circumstellar emission within ~6 AU from Vega is identified at the 3-cr level. The most straightforward scenario 
consists in a compact dust disc producing a thermal emission that is largely dominated by small grains located between 0.1 and 0.3 AU 
from Vega and accounting for 1.23 ± 0.45% of the near-infrared stellar flux for our best-fit model. This flux ratio is shown to vary 
slightly with the geometry of the model used to fit our interferometric data (variations within ±0.19%). 

Conclusions. The presence of hot exozodiacal dust in the vicinity of Vega, initially revealed by K-band CHARA/FLUOR observa- 
tions, is confirmed by our H-band IOTA/IONIC measurements. Whereas the origin of the dust is still uncertain, its presence and the 
possible connection with the outer disc suggest that the Vega system is currently undergoing major dynamical perturbations. 

Key words. Instrumentation: high angular resolution - techniques: interferometric - circumstellar matter 
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^ 1 . Introduction proposed by 'Wilner et al.' ('2002'), and modelled later by IWvattI 

, (2006) and Reche et al. (2008). Infrared observations revealed a 

• The discovery of a debris disc around Vega (HD 172167, AOV, ^^^^^^ axi-symmetric structure extending from 85 to at least to 

> . 7.76 pc) by |Aumannetal.| ((1984) was one of the first hints that gigAU and co ntaining abo ut 3 x lO'^M^ of dust grain s (e.g., 

k> extrasolar planetary systems exist. The debris disc was identi- |suetal.l llool lRieke etal.1 l2005t ISibthorpe et al.1 Mdi) . This 

^ _ fied thanks to an infrared excess (beyond 12^m) with respect ^^^^ ^^^^ ^i^^ ^^^^ j„ j^f^^^^ ^^^^ ^ ^^^^^^ 

^ , to the expected photospheric flux by the Infrared Astronomical ^^^^^^ ^g^g^^j questions about the mechanism at the ori- 



Satellite (IRAS). During subsequent decades, a successive gen- ^-^ ^^^^ IsTet dl (12001 suggest that the extended disc 
eration of instruments (e.g., ISO, SCUBA, Spitzer, AKARI, ^^^^^^ ^ ^^^^^t ^^^^^^^ collision of planetesimals and 
Herschel) confirmed its existence over a wide wavelength range subsequent collisional cascade. This would produce a high- 
covering the infrared and millimetre regimes. The appearance of ^^^^ ^^^^ composed of very small grains (less than a few /.m), 
the disc has been found to vary significantly across this wave- ^^ich are blown out of the system by radiation pressure imme- 
length domain. Sub-miUimetre and milhmetre observations re- ^j^^^jy ^p^^ creation, resulting in the observed large dis c extent, 
vealed a clumpy structure of large dust grains located between Qj^g^ y^g^ ^^^^^^ 455 Mvr. lYoon et aljloTo l). the sta- 
about 80 and 120AU_^g Holland et al. 1998; K oerneretalJ jj^tj^^j nkelihood of such an event occurring with two bodies of 
|2001t IWilneretalJ |2002|; [Marsh et al. , ^006)l!| If present, the ^^^^^^ ^^^^ explain the sub-millimetre observations is , how- 
clumps are thought to be dust grains trapped into mean mo- ^^^^^ ^^^^^ J^itt & D ent 2002). Conversely, MullelitiD 
tion resonances with a planet located near to the disc, as first j20Toh succeeded in reproducing the surface brightness radial 

profile using intermediate size grains in elliptical orbits around 

These results are however challenged by more recent observations ^^^^ planetesimal ring, and therefore conclude that it is 
with the Plateau de Bure Interferometer l lPietu et alJI2011l) . 
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Fig. 1. Array geometry given with the available stations. Stations 
used during the observations are marked with a letter. Telescopes 
A and C can move on stations located along the 35 m north- 
eastern arm, while telescope B can move along the 15 m south- 
eastern arm. 

consistent with a steady-state model. Finally, another interest- 
ing model that was able to reproduce the infrared observations 
suggests t hat the debris disc of Vega is the result of icy planet 
formation dKenvon & BromlevlUo OS). 

These studies of the outer part of the debris disc around 
Vega have been complemented by infrared interferometric ob- 
servations of the inner part. The first attem pt was realised wit h 
the PTI interferometer in dispersed mode (ICiardi et al.ll200ll) . 
The 110-m long baseline and the poor spatial frequency cover- 
age were, however, not well adapted to drawing a clear conclu- 
sion and a simple debris-disc model accounting for 3 to 6% of 
Vega's K-band flux was proposed as the most likely scenario to 
explain the observations. A more recent study using the shortest 
baselines at the CHARA array has derived a K-band flux ratio 
between the stellar photosphere and the debris disc of 1.26 + 
0.28% w ithin the FLUOR field-of-view (about 8 AU in radius, 
lAbsil et al. 2006, 2008). In the N band, the best constraint on the 
thermal emission from warm dust was obtained by nulling inter- 
ferometry, with no resolved emission above 2.1% of the level of 
stellar photosp heric emission at separations larger than 0.8 AU 
(iLiu et al.ll20"0 9). The hot dust grains observed in the inner de- 
bris disc of Vega are believed to derive from collisions between 
larger rocky bodies and/or by the evaporation of comets, as in 
the solar zodiacal disc. The inferred dust populations are, how- 
ever, much hotter, more massive, and composed of much smaller 
grains than the zodiacal cloud. Such grains would be expected 
to be expelled from the inner planetary system by radiation pres- 
sure within only a few years, which indicates inordinate replen- 
ishment rates. 

The present paper reports on new infrared interferometric 
observations of the hot inner part of the debris disc around 
Vega. Visibility and closure phase H-band measurements were 
obtained in June 2006 with the 3-telescope IOTA/IONIC in- 
strument, one year after the pioneering K-band observations of 
[Absil et al. (2006). Studying the inner debris disc around Vega is 
of prime importance as it corresponds to the location where ter- 
restrial planets are supposed to be formed and evolve, and will 
provide additional constraints for discussing the still-debated 
origin of the dust. On a longer time scale, the characterisation of 
inner debris discs is also relevant for preparing the programme 
of future space missions dedicated to the direct detection and 
characterisation of Earth-like planets, since the presence of large 
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Fig. 2. Sampling of the Fourier (u,v) plane for the selected data 
set. 



quantities of warm dust in the habitable zone around nearby 
main sequence stars migh t jeopardize the succ ess of such mis- 
sions dRoberge et alj|2009l: iDefrere et al.ll2010h . 

2. Observations and data reduction 

2.1. Instrumental setup 

The interferometric data discussed in this paper were ob- 
tained with the IONIC combiner at the IOTA int erferometer 
(Infrared-Optical Telescope Arrav. lTYaub et alj2003h . IOTA was 
a three telescope interferometer located at the Fred Whipple 
Observatory atop Mount Hopkins (Arizona, USA). The IOTA 
telescopes were movable among several stations along an L- 
shaped track: telescopes A and C can move along the 35-m 
north-eastern arm, while telescope B moves along the 15-m 
south-eastern arm, see Figure [T] This enables synthesizing an 
aperture of 35 m x 15 m, which corresponds to an angular reso- 
lution of about 5-12mas at 1.65yum. The collecting optics con- 
sisted of 0.45-m Cassegrain primary mirrors, which were fed by 
siderostats. Tip-tilt servo systems mounted behind the telescopes 
compensated for the atmospherically induced motion of the im- 
ages, while delay lines actively tracked the fringes by adjusting 
the optical path delay between the different baselines. All the 
beams were deflected by a series of mirrors into the laboratory, 
where they were coupled into single-mode fibres. The fibres then 
fed the spatially filtered beams into the IONIC3 integrated optics 
beam combiner dBereer et alJl2003h . which combined the beams 
coaxially and pairwise with a ratio of 50:50. For each baseline, 
the beam combination produced two complementary outputs, 
which were shifted in phase by n with respect to each other 
and which were recorded on a PICNIC camera ( Pedretti et al] 
2004). Although the information recorded by these two channels 
is in principle redundant, it is used to remove residual photomet- 
ric fluctuations simply by subtracting the signals from the two 
channels. The PICNIC detector array is a 256 x 256 pixel ar- 
ray, arranged in four quadrants of 128 x 128 pixels and sensitive 
to the 0.8-2.5 fim wavelength region. The pixel readouts were 
performed in a non-destructive way. 
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Interferometric observations of Vega were obtained in the 
infrared H-band (1.50 - 1.80//m) with IONIC, using the SE05- 
NE07-FIX and NE10-SE15-FIX triplets of the IOTA array 
(see Figure [1] and the corresponding u-v plane in Figure |2|. 
Observations took place in June 2006, on June 9, June 10 
(NE07-SE05-FIX triplet), June 13 and June 14 (NE10-SE15- 
FIX triplet). A fringe-tracking algorithm was applied in real time 
to ensure that interference occ urs nearly simultaneously for all 
baselines ( lPedrettietalJl2005h . In addition, the light was split 
into two orthogonal polarisation axes with a Wollaston prism af- 
ter beam combination, and these were recorded simultaneously. 
This choice was initially made to improve the stability of the 
instrumental response, hence the data quality. However, a non- 
linear regime of the detector, probably related to partial satura- 
tion effects, was observed for several pixels in the data from June 
9 and June 10. This non-linearity arose only on the pixels cor- 
responding to the first polarisation, which presents significantly 
higher fluxes than the second one. Consequently, the first polar- 
isation in the data from June 9 and June 10 has been discarded. 
This behaviour has not been observed in the data of June 13 and 
June 14 in any of the two polarisations. 

2.2. Data reduction 

The raw squared visibilities and closure phases w ere obtained by 
using the es tablished IDE routines described by iMonnier et al.l 
(I2004l2006b . The interferometric transfer function of the instru- 
ment, i.e., the response of the system to a point source, is esti- 
mated with calibrator stars. All calibrator stars (listed in Table|2]l 
were ch osen from tw o catalogues developed for this specific 
purpose dBorde et al.ll2002: Merand et al..,2005.) . The statistical 
and systematic error bars on the transf er function est i matio n 
were computed following Appendix C of iKervella et al.l (l2004l) . 
Based on the calibrator measurements, the interferometric trans- 
fer function was e stimated for the whole night by polynom ial in- 
terpolation (as in Lebouquin et alJl2006tlAbsil et al.ll2009 ). The 
raw squared coherence factors ff for Vega and its calibrators 
are shown in Figure[3]for a representative night. A second-order 
polynomial was found appropriate to following the variations in 
the interferometric transfer function during the whole night as 
the target and calibrator moved across the sky. The statistical er- 
ror bar on the interferometric transfer function was computed 
locally at any given time by using the covariance matrix on the 
parameters of the polynomial fit. The systematic error bar, on the 
other hand, was computed globally with a weighted sum of the 
systematic error bars on all transfer function estimations, taking 
the correlation between calibrators into account. 

Since all calibrators chosen in this study are late K giants, 
the squared visibilities and closure phases were computed for 
the wide spectral bandwidth on which the IOTA observations 
were performed, considering the actual spectrum of the star (us- 
ing tabulated H-band spectra from Pickles 1998) and the spec- 
tral transmission of the IOTA/IONIC instrument. An additional 
systematic error was considered in order to consider the chro- 
maticity effects that limit the absolute precision of calibration. 
For squared visibilities, we estimated the impact of chromatism 
using spectrally dis persed data that were ob tained the same year 
with IONIC (see iLacour et alJ 120081: iPedretti et all 12009 ). Our 
analysis shows that the transfer function does not vary signifi- 
cantly with the spectral type of the observed star for any given 
instrumental setup, but that the shape of the transfer function 
with respect to wavelength depends on the instrumental set-up 
(baseline, detector, read-out mode, etc.) in an unpredictable way. 
Therefore, we cannot accurately predict the colour correction 
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Fig. 3. Squared coherence factors (/i^) and interferometric trans- 
fer function (TF) estimations (T^) for one representative observ- 
ing night (June 9th) on the NE07-SE05 baseline. Filled symbols 
correspond to //^ data points for Vega (red diamonds) and the 
calibrator targets (blue squares). The interferometric TF is rep- 
resented by the empty blue squares, which are derived from the 
calibrator measurements taking their known diameter into ac- 
count. The solid line is an interpolation of the TF, obtained by 
fitting a second-degree polynomial to its estimated values. The 
1-cr statistical and systematic error bars on the estimation of the 
TF are, respectively, represented by dotted and dashed lines. The 
two gaps during the night (one around 8:00 UT and the other 
around 9:30 UT) are due to the change in the long delay line 
offset and to some passing clouds. 

factor for V^ measurements between the various observed spec- 
tral types. This variable behaviour of the spectral transfer func- 
tion across the H band has an impact on the calibration accuracy 
of + 1 % at most between an AOV and a K3III star, which we add 
quadratically to our estimated error bars. For closure phases, we 
used a conservative systematic error bar of 0.5° to account for 
chromatism effects. This value is based on engineering tests car- 
ried out by (Monnier et al. 2 006) . which showed that the closure 
phase presents a systematic error of this order between a hot star 
(B8) and a cool star (M3). 

Finally, for the nights of June 13 and 14 (large triplet), 
only data corresponding to hour angles ranging between -2 and 
2 hours have been used in the following discussions to achieve 
good accuracy on the calibrated squared visibilities (see indi- 
viduals measurements in Table [T]i. The final calibrated data are 
presented in Figure |4] (squared visibility) and Figure |5] (closure 
phase). 

3. Modelling the data with simple models 

As shown bv lAbsilet all (|2006|) wifli CHARA/FLUOR, a real- 
istic stellar photospheric model cannot reproduce K-band vis- 
ibility measurements. As shown in Figure |4] this is also the 
case for our H-band IOTA/IONIC data with measured visibili- 
ties clearly below the expected level of the photospheric visibil- 
ity. These observations cannot be reconciled with a purely pho- 
tospheric model, because on the one hand short-baseline data 
are weakly sensitive to the model parameters (angular diameter 
and limb-darkening profile), and on the other, the mo del param- 
eters are already known w ith good accuracy ( Aufd enberg et al.l 
120061: iPeterson et al.ll2006l) . Furthermore, the visibility deficit is 
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Table 1. Individual measurements obtained with IOTA/IONIC in June 2006. Only one polarisation is kept in the data of June 9 and 
10, while both polarisations are kept for June 13 and 14. 





UT 


HA 


# files 


Array 




Proj. baseline [m] 






Date 


[lih:mm] 


[hh:mm] 


rec. 


pre. 


configuration 


AB 


BC 


AC 


CALi 


CALz 


2006/06/09 


5:47 


-3:04 


2 


2 


NE07 SE05 FIX 


8 


22 


4.83 


5.86 


0Her 


fHer 




6:03 


-2:48 






NE07 SE05 FIX 


8 


29 


4.86 


6.04 


0Her 


eHer 




6:12 


-2:38 






NE07 SE05 FIX 


8 


32 


4.87 


6.14 


SHer 


n Her 




6:29 


-2:22 






NE07 SE05 FIX 


8 


37 


4.89 


6.30 


n Her 


^Lyr 




6:44 


-2:07 






NE07 SE05 FIX 


8 


41 


4.92 


6.43 


iLyr 


^Lyr 




7:02 


-1:49 






NE07 SE05 FIX 


8 


45 


4.94 


6.57 


^Lyr 


/iLyr 




7:17 


-1:34 






NE07 SE05 FIX 


8 


48 


4.96 


6.67 


^Lyr 


/iLyr 




7:28 


-1:23 






NE07 SE05 FIX 


8 


50 


4.98 


6.73 


Ahyr 


^Lyr 




7:39 


-1:12 






NE07 SE05 FIX 


8 


51 


4.99 


6.79 


iLyr 


^Lyr 




8:48 


0:03 






NE07 SE05 FIX 


8 


56 


4.99 


7.01 


n Her 


n Her 




8:58 


0:07 






NE07 SE05 FIX 


8 


57 


4.98 


7.02 


n Her 


n Her 




9:06 


0:16 






NE07 SE05 FIX 


8 


57 


4.97 


7.03 


n Her 


n Her 




9:15 


0:25 






NE07 SE05 FIX 


8 


57 


4.95 


7.04 


71 Her 


0Her 




9:24 


0:34 






NE07 SE05 FIX 


8 


57 


4.93 


7.04 


SHer 


0Her 




10:30 


1:40 






NE07 SE05 FIX 


8 


56 


4.68 


7.01 


0Her 


K Lyr 




10:40 


1:50 






NE07 SE05 FIX 


8 


55 


4.62 


7.00 


K Lyr 


/iLyr 




10:57 


2:07 






NE07 SE05 FIX 


8 


54 


4.51 


6.98 


K Lyr 


iLyr 




11:06 


2:16 






NE07 SE05 FIX 


8 


53 


4.45 


6.97 


K Lyr 


iLyr 




11:17 


2:27 






NE07 SE05 FIX 


8 


52 


4.37 


6.96 


^Lyr 


/iLyr 




11:26 


2:36 






NE07 SE05 FIX 


8 


51 


4.30 


6.95 


^Lyr 


K Lyr 




11:37 


2:47 






NE07 SE05 FIX 


8 


50 


4.21 


6.94 


K Lyr 


K Lyr 


2006/06/10 


8:10 


0:37 


3 


2 


NE07 SE05 FIX 


8 


55 


5.01 


6.93 


SHer 


0Her 




8:20 


0:27 


3 


3 


NE07 SE05 FIX 


8 


55 


5.00 


6.96 


0Her 


eHer 




8:31 


0:16 


3 


3 


NE07 SE05 FIX 


8 


56 


5.00 


6.98 


0Her 


eHer 




8:42 


0:05 


3 


3 


NE07 SE05 FIX 


8 


56 


4.99 


7.00 


SHer 


0Her 




8:52 


0:06 


3 


3 


NE07 SE05 FIX 


8 


57 


4.98 


7.02 


SHer 


0Her 




9:06 


0:19 


3 


3 


NE07 SE05 FIX 


8 


57 


4.96 


7.03 


0Her 


eHer 




9:16 


0:30 


3 


2 


NE07 SE05 FIX 


8 


57 


4.94 


7.04 


0Her 


n Her 




9:32 


0:46 


3 


3 


NE07 SE05 FIX 


8 


57 


4.90 


7.04 


71 Her 


7t Her 




9:45 


0:58 


3 


3 


NE07 SE05 FIX 


8 


57 


4.86 


7.04 


71 Her 


7t Her 




9:55 


1:09 


3 


3 


NE07 SE05 FIX 


8 


57 


4.82 


7.03 


Tc Her 


n Her 


2006/06/13 


9:35 


1:01 


2 


2 


NEIO SE15 FIX 


17 


72 


14.56 


10.01 


0Her 


K Lyr 




9:47 


1:13 


2 


1 


NEIO SE15 FIX 


17 


66 


14.42 


10.00 


K Lyr 


K Lyr 




10:01 


1:27 


2 


1 


NEIO SE15 FIX 


17 


58 


14.24 


9.98 


K Lyr 


0Her 




10:28 


1:53 


2 


2 


NEIO SE15 FIX 


17 


40 


13.81 


9.95 


K Lyr 


eHer 


2006/06/14 


7:08 


-1:23 


2 


2 


NEIO SE15 FIX 


17 


99 


14.92 


9.57 


7X Lyr 


n Her 




7:19 


-1:12 


2 


2 


NEIO SE15 FIX 


17 


98 


14.96 


9.65 


71 Lyr 


71 Her 




7:31 


-0:60 


2 


2 


NEIO SE15 FIX 


17 


98 


14.98 


9.73 


71 Lyr 


0Her 




7:51 


-0:40 


2 


2 


NEIO SE15 FIX 


17 


96 


15.01 


9.84 


0Her 


eHer 




8:01 


-0:29 


2 


2 


NEIO SE15 FIX 


17 


95 


15.01 


9.88 


SHer 


K Lyr 




10:23 


1:52 


2 


2 


NEIO SE15 FIX 


17 


41 


13.83 


9.95 


SHer 


0Her 



Cols. 1, 2, and 3 give the date, UT time, and the hour angle, cols. 4 and 5 give the corresponding number of files recorded and kept during the data 
reduction process, col. 6, 7, 8, and 9 give the array configuration and corresponding projected baselines, respectively. The last two columns give 
the calibrator stars used before (CALi) and after (CAL2) the observation of Vega. 

Table 2. Fundamental parameters and estimated angular diameter for Vega {a Lyr) and its caUbrators. 



Identifier HD number RA-J2000 DEC-J2000 Sp. type my m„ Gw ± Icr Refs. 

[d m s] [dm s] [mas] 



Ti Her 


156283 


17 15 02.83 


+36 48 32.98 


K3II 


3.16 


0.12 




[1,2] 


0Her 


163770 


17 56 15.18 


+37 15 01.94 


Kllla 


3.85 


1.14 


3.15»"3'' 


[1,2] 


K Lyr 


168775 


18 19 51.71 


+36 03 52.37 


K2IIIab 


4.32 


1.86 


2 280-025 


[1,2] 


a Lyr 


172167 


18 36 56.33 


+38 47 01.29 


AOV 


0.03 


-0.03 


3 3120.067 


[1,3] 


iLyr 


176670 


19 00 00.83 


+32 08 43.85 


K2.5III 


4.95 


1.84 


2 410.026 


[1,2] 



Refere nces. Co ordinates, spectral types and magnitudes from [1] SIMBAD; limb-darkened diameters and 1-cr errors from [2] |Bordeetai] ( l2002l) . 
and [3] lAbsiletal. ( 2008) . 

observed on the three different baselines with roughly the same ness of Vega's photosphere (which is not significant since Vega 
magnitude so that it cannot be explained by the apparent oblate- is seen almost pole-on). 
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To explain this visibility deficit. lAbsil et alJ (l2006l) proposed 
a model of a star surrounded by an exozodiacal disc accounting 
for 1.26 + 0.28% of the K-band stellar flux within the FLUOR 
field-of-view (~8AU in radius. IXbsil et all I2b08h . In particu- 
lar, they addressed and ruled out a series of scenarios. For in- 
stance, a bound low-mass companion was also shown to be a 
viable explanation although much less probable. Other expla- 
nations, such as a bright background object within the FLUOR 
field-of-view (probability less than 10"^) or stellar spots (appear- 
ing in the second and higher lobes of the visibility function), 
were ruled out with good confidence. At that time, it was argued 
that stellar winds and mass-loss events could also be the origin 
of the near-infrared exce ss, but thes e scena rios were shown to 
be very unlikely later bv lAbsil et al.l (l2008l) . In the two follow- 
ing sections, we consider the two most probable interpretations 
separately (the exozodiacal disc and the binary companion) and 
check whether simple models can reproduce our measurements. 

3.1. The exozodiacal disc scenario 

To constrain the near-infrared excess emission around Vega, 
three simple models of a circumstellar disc have been fitted to the 
data: a diffuse emission uniformly distributed across the field- 
of-view, the zodiacal disc model of Kelsall et al. ( 1998), which 
is implemented in the Zodipic package^ and a narrow ring of 
dust located at twice the sublimation radius (rsub=0.1 AU for a 
sublimation temperature Tsub=1700K). All models are assumed 
to be point-symmetric, as suggested by the calibrated closure 
phases (see Figure|5j, and use an effective temperature of 8027 K 
for Vega as viewed by the dust in the equatorial plane (Jason 
Aufdenberg, private communication). The results of this fitting 
procedure are shown in FigureH) with the corresponding H-band 
disc/star flux ratios and reduced in Table [3] All three models 
fit the data set (xj. =0.5-0.7) equally well with a best-fit flux ratio 
of 1.35 + 0.49%, 1.26 + 0.45%, and 1.20 + 0.43%, respectively, 
for the uniform disc model, the zodiacal disc model, and the ring 
model. All the best-fit contrasts are compatible within 1 cr, which 
indicates that the final flux ratio does not heavily depend on the 
distribution of the circumstellar emission. 

To confirm this conclusion, the squared visibilities from syn- 
thetic images of a geometrically and optically thin debris disc 
were computed and compared to the data set. The synthetic im- 
ages are based on three parameters: the disc/star flux ratio, the 
inner disc radius rin, and the exponent a of the power law de- 
scribing the density decrease as a function of distance (rir oc r"). 
The synthetic images only include thermal emission, which is 
expected to be largely dominant for the dust temperatures ex- 
plored here, and they assume pure blackbody emission (grain 
temperature proportional to r""^, with T = 1700 K at 0.1 AU). 
For each couple of parameters (ri„,a) ranging between the sub- 
limation radius and 2 AU and between and -50, respectively, 
the best-fit flux ratio was computed by optimisation. The re- 
sults show that all models fit the data equally well with a reduced 

ranging between 0.48 and 0.84, confirming exozodiacal dust 
as a good scenario to explain the deficit of visibility. It is, how- 
ever, impossible to conclude on the location of the dust based on 
the sole interferometric data. The best-fit disc/star flux ratios are 
all compatible within approximately 1-cr error bars and range 
from 1.04% to 1.42%. We therefore adopt a H-band flux ratio 
of 1.23 + 0.53% in the following discussions. The error bar in- 



Vega (H bard, v.urie 2006) 



^ Zodipic is an IDL program developed by M. Kuchner and C. Stark 
for synthesizing images of exozodiacal clouds. It can be downloaded 
from http://asd.gsfc.nasa.gov/Marc.Kuchner/home.htmli 
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Fig. 4. Expected squared visibility as a function of the baseline 
length for Vega's limb darkened photosphere (blue solid line), 
given with the measured squared visibilities and related error 
bars (blue symbols for June 9 and 10, green symbols for June 

13 and 14 on polarisation 1, and red symbols for June 13 and 

14 on polarisation 2). The thickness of the blue solid line cor- 
responds to the 5-cr uncertainty related to the uncertainty on the 
stellar diameter The best-fitted models of a limb-darkened pho- 
tosphere surrounded by a uniform circumstellar emission (dot- 
ted blue line), a zodiacal disc (dashed green line), and a narrow 
ring (dash-dotted red line) are also represented (see discussion 
in Section 3). Residuals of the fit are given in the bottom panel 
for the uniform disc model. 

Table 3. Best-fit flux ratio and goodness of fit for 3 different disc 
models. 



Disc model 


Flux ratio [%] 


xi 


Uniform emission 


1.35 


±0.49 


0.52 


Zodiacal disc 


1.26 


± 0.45 


0.52 


Narrow ring [0.2 AU] 


1.20 


±0.43 


0.62 



The angular diameter and effective temperature of the stellar photo- 
sphere have been fixed respectively to 3.312 ± 0.067 mas (see Table|2) 
and 8027 K as viewed by the dust in the equatorial plane (see main text). 



eludes both the uncertainty due to the model (+0.19%) and the 
maximum error bar on an individual model (+0.49%). This in- 
ferred disc/s tar contrast is marginally compatible with the one 
predicted bv lAbsil et all 12006) based on the CHARA/FLUOR 
observations (about 0.6% with an inner rim located between 0.17 
and 0.3 AU). The implications on the disc properties of this new 
H-band measurement are discussed in Section|4| 

3.2. Tiie binary scenario 

The possible presence of a low-mass binary companion around 
Vega has been the subject of various studies (e.g., radial veloci- 
ties, astrometric measurements, and high dynamic range single- 
pupil imaging), drastically reducing the range of allowed param- 
eters for the putative companion. Before discussing the existing 
data in the literature, it is necessary to check whether a binary 
companion could effectively reproduce our interferometric ob- 
servations. Given the relatively good sampling of the u-v plane 
(see Figure |2]i and the uniform deficit measured along the dif- 
ferent baselines (see Figure|4]i, we might expect to derive strong 
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Fig. 5. Calibrated closure phase measurements given with error 
bars as a function of hour angle. The color code is the same 
as in Figure |4] The mean closure phase is represented by the 
black solid line, while the dotted lines represent the 1-cr statisti- 
cal uncertainty on its value. The bottom plot shows the residual 
between the data and the mean value of closure phases. 



constraints on the properties of a binary companion based on the 
sole interferometric measurements. 

To assess whether a low-mass binary companion would be 
compatible with our observations, we focus as a first step on 
the calibrated closure phases, which are shown in Figure |5] as a 
function of the datafile number The mean closure phase is 0.74°, 
with a 1-cr error of +0.73°. This error includes the statistical er- 
ror due to the scatter in the individual CP measurements and 
the systematic error due to the chromatism of the beam com- 
biner (0.5° as discussed previously), which clearly contributes 
to the low positive value of the mean closure phase. Basically 
all closure phases are consistent with the mean value (x^. of 
0.69) for all the various projected baselines of both configura- 
tions, suggesting that any deviation from point symmetry in the 
near-infrared excess source is not larger than 1 % for the spatial 
frequencies explored here. To determine which combination of 
orbital parameters and contrast of the binary companion would 
be compatible with these measurements, we fitted a model of a 
binary star to the whole data set (squared visibilities and clo- 
sure phases) in a second step. The orbital motion of the puta- 
tive stellar companion is ta ken into accou nt assuming a mass of 
2.3 Mq for Vega ( Aufdenber g et al.ll2006l) . The mass of the com- 
panion is com puted using the evolutionary models developed by 
iBaraffe et al.l ( lI998), assuming it h as the same age as Vega itself 
(about 455Mvn lYoon et al.ll2010l) . We also assumed the com- 
panion to be orbiting on a circular orbit within the plane of the 
sky (Vega is almost se en pole-on with an inclination angle esti- 
mated to 5.7°. lAufdenberg et al.ll2006 ). The model is therefore 
based on three parameters: the semi-major axis, the orbital phase 
at a given time to, and the binary flux ratio. For each couple of 
semi-major axis/contrast, we computed the orbital phase at to 
that minimizes the xj- With this, we have produced x^- maps, 
which are represented in Figure |6] as a function of semi-major 
axis and binary contrast for the squared visibilites (left) and the 
closure phases (right). 

The x^ maps show that a range of models fits the interfer- 
ometric data in a satisfactory way, the effect of the Gaussian 
beam profile being clearly evident in the lefthand figure as the 
binary contrast required to fit the data increases with the bi- 



nary separation. All models corresponding to a binary separa- 
tion larger than 1.5 AU can be rejected with a 5-cr confidence 
level by the closure phase map. Within 1.5 AU, the squared visi- 
bility map shows that only models corresponding to a contrast 
between 0.6 and 1.6% are compatible with the data with the 
same confidence level. Furthermore, a very short semi-major 
axis (<0.1 AU) can be rejected because it is either incompati- 
ble with the squared visibilities or with the closure phases, re- 
gardless of the contrast. Therefore, based on our interferometric 
data, the only possible parameters for a putative binary compan- 
ion are a semi-major axis ranging between 0.1 and 1.5 AU and 
a contrast ranging between 0.6 and 1.6%. Such properties for 
the companion are, however, rejected by Hipparcos astrometric 
measurements of Vega, which have shown that the orbital semi- 
major axis of a putative companion cannot be larger than 6.3 mas 
(=0 .05 AU=4R^^) with a 99% confidence ass uming a circular or- 
bit (iPerrvman etal]|1997l: IXbsil et all 12006 ). The absence of a 
detected companion in radial velocity measurements confirms 
this conclusion further. All these elements suggest that the pres- 
ence of a close companion around Vega to explain the observed 
visibility deficit is extremely unlikely. 

4. Further constraints on the exozodiacal disc 

The analysis presented in the previous section confirms that the 
most straightforward scenario for explaining the visibility deficit 
is the presence of hot dust grains producing thermal emission 
in the near-infrared. To constrain the disc parameters, several 
archival spectro-photometric and interferometric measurements 
at near- and mid-infrared wavelengths exist in the lit erature (see 
Table jH ). In addition to the measurements used bv lAbsil et al.l 
(I2006h . we add our direct measurement of the H-band excess 
flux (1.23 + 0.53%) and the IRS spectrum found in the Spitzer 
archives. The Spitzer spectrum has been reduced with the c2d 
pipeline developed by iLahuisI (12006) . and binned into a few 
equivalent broad-band photometric measurements for the sake 
of spectral energy distribution (SED) modelling. The idea is then 
to check whether there is at least a disc model that is compatible 
with all the constraints and derive its properties. 

To reproduce the SED of the infrared excess as listed in 
Table |4] radiative trans fer modelling was performed using the 
code developed by Au gereau et al.l (Il999h for cold d ebris discs, 
and adapted to the case of exozodiac al discs (e.g.. lAbsil et al.l 
2006; di Folco et al. 2007t lAbsil et al.. .2008. 2009). The code 
considers a population of dust grains with a parametric surface 
density profile and size distribution, and computes the thermal 
equilibrium temperature of the grains exposed to the stellar ra- 
diation. It can handle various chemical compositions for the 
grains in order to investigate for instance the fraction of carbona- 
ceous or silicate material contributing to the observed emission. 
Specific care is given in the model to the treatment of exozodia- 
cal dust close to the sublimation radius, to account for the size- 
dependent position of this radius. On output, the code computes 
the thermal and scattered light emissions of debris disc mod- 
els and produces both images and SEDs over a broad param- 
eter space to combine observations from different instruments 
(IOTA/IONIC, CHARA/FLUOR, MMT/BLINC, Spitzer/IRS, 
broad-band fluxes), taking their specific response into account. 

Using the code described above with an effective tempera- 
ture for Vega of 8027 K as viewed by the dust in the equatorial 
plane, a grid of debris disc models was computed with the fol- 
lowing parameters: the peak position ro of the surface density 
profile (possible values ranging from 0.05 to 1 AU, provided that 
the sublimation radius is not reached), the minimum size Omin of 
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Fig. 6. Reduced map for the fit of a simple binary model to the squared visibility (left) and closure phase (right). The maps are 
given as a function of two free parameters: the semi-major axis of the companion's orbit and the flux ratio between the putative 
companion and Vega. For each couple of semi-major axis/contrast, the orbital phase at a given time to has been optimised to 
minimize the x^- All models corresponding to a semi-major axis larger than 2 AU in the closure phase map are beyond 5 cr and are 
not shown for the sake of clarity. 



Table 4. Available constraints on the near- and mid-infrared ex- 
cess around Vega. 



Wavelength 


F 

A meas 


Excess 


Instruments 


Ijim] 


[Jy] 


[%] 




1.26 


1574 + 34 


2.4 ± 2.9% 


Catalina',UKIRT" 


1.60 


1055 + 32 


-2.4 ±3.6% 


Catalina' 


1.65 




1.23 ±0.53% 


IOTA/IONIC 


2.12 




1.29 ±0.19% 


CHARA/FLUOR' 


2.12 




5+i% 


p'pj4 


2.20 


655 ± 14 


4.1 ±"3.0% 


Catalina', UKIRT^ 


3.54 


283 ± 6 


3.1 ± 3.0% 


Catalina', UKIRT^ 


4.80 


167 ± 8 


7.1 ± 5.1% 


Catalina', UKIRT- 


10 


40.0 ± 3 


6.0 ± 4.5% 


Various^ 


10.5 


33.7 ± 0.33 


0.0 ± 5.0% 


IRS spectrum 


10.6 




0.2 ± 0.7% 


MMT/BLINC'' 


12 


27.0 ± 0.3 


1.2 ± 5.0% 


IRAS^ 


12.5 


23.89 ± 0.39 


0.0 ± 5.0% 


IRS spectrum 



The superscript on the name of the ins trum ent in the last co l umn g ives 
the reference: (l)'CamDins et al.' ("1985), ('2i'Blackwell et alJ (ll983h . (3) 
Absil et al. (2008), (4)" Ciardi et al. (2001), (5) Rieke et alj ([1985^, (6) 
Liu etalJ (120041 (7),Cohen et al. (1992), with the absolute photometric 
error estimated by|Aumann et al. ( 1984). The interferometric data from 
IONIC, FLUOR, PTI, and BLINC only sample a specific part of the 
inner disc, while the photometric studies include Vega's entire environ- 
ment. 



the size distribution (from 0.01 to 54.6/vm), the slope k of the 
size distribution (from -2.7 to -9), the slope a of the surface den- 
sity profile beyond ro (from to -9.5), and the volume fraction C,- 
of carbon grains (from to 100%). Each model assumes no az- 
imuthal dependence and sublim ation temper atures of 1 200 K for 
silicate grains ("astrosilicates", Drainel2p03[) and 1900 K fo r car- 
bonaceous grains ("amorphous carbon". IZubko et"al]|1996h . The 
size distribution (between a^in and flmax) has been accordingly 
truncated at the sublimation size (asub), which depends on the 
radial distance to the star since only grains large enough to sur- 
vive the sublimation process can actually coexist (see Figure|7]i. 



For each model on the grid, the SED has been computed 
and compared to the measurements reported in Table |4] The 
first obvious result obtained from the xj- analysis is that pure 
silicate models (Cr-Q) do not fit the SED well with a x^r ^1" 
ways larger than 7. The goodness of the fit improves signif- 
icantly by introducing carbonaceous grains, which can sur- 
vive closer to Vega due to their higher sublimation tempera- 
ture. This conclusion is supported by previous studies explain- 
ing the lack of significant silicate emission features around 
lOfim by the presence of large amounts of highly refractive 
grains in the inner disc (Laor & Drainel[T993t IZubko et al.lll996t 
[Gaidos & Koresko 2004). The actual amount of carbonaceous 
grains can, however, not be constrained since the xj- typically 
varies from 1.3 for a carbonaceous-poor disc (Cf=5%) to 0.8 
for a pure carbonaceous disc (C, = 100%) in the steep power law 
regime (a < -3). For a flatter power law {a > -3) like in our zo- 
diacal disc (a - -0.34, Kelsall et al. 1998), the disc models do 
not fit the various measurements well with a x^ always higher 
than 5 even for pure carbonaceous discs. In addition, relatively 
small grains (omin < 1 jum) with a steep size distribution profile 
(k <-3) are also required to obtain a good fit to the SED (x^ < 3). 
Two good-fit models are represented in Figure|2]by flux density 
maps given as a function of grain size and distance to the star. 
The upper figure corresponds to the best-fit model (x^ = 0.84, to- 
tal mass = 1 .9 X 10"'' M©) obtained for Omin - 0.2 fim and k = -5, 
while the bottom figure corresponds to the same model with Oniin 
equal to 0.01 jum Or? = 1.01, total mass = 2.1 x lO^'^M®). The de- 
pendence of sublimation temperature on the grain size is clearly 
visible in these figures, the larger grains surviving closer to Vega 
than the small ones. The figures also show that the largest grains 
located near 0. 1 AU contribute relatively less to the total flux 
than the smaller grains located further away from the star and 
present in larger numbers. In both cases, the flux density de- 
creases very quickly with the distance to Vega as shown by the 
various contours. This geometry is mostly constrained by the in- 
terferometri c measurements , and particularly the non-detection 
reported bv lLiu etal.! (l2004l) with the MMT/BLINC, which re- 
quires a steep power law in order to reduce the amount of dust in 
the regions further than 1 AU. In addition, recent observat ions of 
Vega with the Palomar Fiber NuUer (~2.2 (im. .Mennesson et al.l 
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Fig. 7. Flux density maps given as a function of grain size and 
distance to the star computed for Oniin = 0.2 yum (best-fit model, 
upper figure) and a^nm = 0.01 jum (lower figure). A surface den- 
sity power law of a = -3 and a size distribution power-law ex- 
ponent of K - -5 have been assumed in these plots, with a maxi- 
mum size of 1000 yum. In this simulation, the disc is composed of 
50% silicates and 50% carbonaceous grains (C,.=50%, see text). 
The contours are plotted every power of 10 between 10 ' and 
10"^ of the maximum flux density. 



suggest that any dust population contributing to at least 
1% of the near-infrared excess can arise only within 0.2 AU, a 
conclusion in favour of a^nin equal to 0.2yum in order to suffi- 
ciently reduce the dust emission beyond 0.2 AU. Finally, we used 
a Bayesian approach to compute the normalised probability den- 
sity of the disc models based on the grid. The results of this 
analysis shows that the disc mass has a probability density that 
peaks at about lO^'^M(g (see Figure |9]), equivalent to the mass of 
an asteroid about 20 km in diameter (assuming 2.5 g/cnr' den- 
sity). 

In summary, the most straightforward scenario consists in a 
compact debris disc with a steep density profile (a < -3) pro- 
ducing a thermal emission that is largely dominated by small 
(flmin < l.Oyum) silicates and carbonaceous grains located be- 
tween 0. 1 and 0.3 AU from Vega. A representative SED is shown 
in Figure [8] together with the SED of the stellar photosphere 
as seen pole-on (upper solid curve, Teff= 10150 K) and as seen 
by the exozodiacal dust (lower solid curve, Tett=8027K). This 
model assumes a size distribution exponent of -5 with limiting 
grain sizes of a^i^ - 0.2//m and a^ax - 1000 yum, a surface den- 
sity power law of -3, and a disc composed of 50% silicates and 
50% carbonaceous grains. By using this disc model to fit our in- 
terferometric data set, we derived a best-fit contrast of 1.23 + 
0.45%. Finally, the flux ratio derived here is about 1.5 o" above 
the value extrapolated from the K-band CH ARA/FLUOR mea - 
surements (expected H-band excess of 0.6%. lAbsiletani2006h . 




10"'° 10"^ 10"" 10"^ 10"*^ 

M„„„ [M,„„] 



Fig. 9. Normalised probability density for the disc total mass 
computed with a Bayesian approach based on the grid. The 
probability density peaks at about 10"'' M®, equivalent to the 
mass of an asteroid about 20 km in diameter. 



Although this model successfully reproduces the Vega SED, in- 
cluding both CHARA/FLUOR and IOTA/IONIC flux ratios, an- 
other possible scenario that could explain the small discrepancy 
(at least partially) is that the brightness of the exozodiacal disc 
around Vega is evolving on time scales as short as one year. This 
issue is discussed briefly in the next section and is currently the 
subject of follow-up oservations of Vega at the CHARA array. 

5. Discussion 

Several scenarios can explain the presence of hot dust in 
the close vicinit y of Vega. The age of Vega (about 455 Myr, 
lYoon et al]|2010h precludes the possibility for this dust to be pri- 
mordial, because the time scale to remove suc h dust is shorter 
than 10 Myr (Bac kman & Parescel[T99l IWvatt 2008). The de- 
tected dust is therefore second generation and produced either 
by the collision of larger bodies or by cometary sublimation. 
In addition, a large dust production rate (~ 10"^Me/yr) is nec- 
essary to account for the amount of detected dust, which has 
a limited lifetime (of a few years at most) in the inner Vega's 
system due t o radiation pressure and colhsional destru ction 
(iKrivovetal .1 12000). As pointed out by lAbsil et al] ( l2006l) . the 
inward transport of dust from the outer disc (located at about 
85 AU, e.g., iSibthorpe et al. 201 0) through P-R drag is very in- 
efficient given t he long time scale of this process (2x10^ yr, 
iDent et al]l2000h and the relatively short collisional time scale 
in the outer disc (5x10^ yr). Instead, the dust would be gener- 
ated locally by comets orbiting in the inner system. The comets 
would likely originate from the outer disc of Vega or possibly 
from an inner population of icy bodies as in the case of /3 Pic 
(iBeust & Morbidelli''2000'). In fact, comets are believed to be the 
origin of at least 90% of the dust observed in the solar zodiacal 
cloud dNesvornv et al.l 120101) . Because they are exposed to the 
stellar radiation and wind, the dust particles released in the inner 
Vega system could then pile up in a zo ne close the sublimatio n 
radius with a very steep density profile (iKobavashi et al.ll2009l) . 

In any case, the high production rate needed to account for 
the amount of detected dust suggests that Vega is currently un- 
dergoing major dynamical perturbations. A dynamical "shake- 
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Wavelength, X Cjxm) 

Fig. 8. A possible good fit of our debris disc model to the photometric and interferometric constraints of Table |4] (y^ - 0.84). The 
model used here has a size distribution exponent of -5 with limiting grain sizes of amin = 0.2 fim and a^ax = 1000 yum, a surface 
density power law of -3, and it assumes a disc composed of 50% silicates and 50% carbonaceous grains. The dashed and dotted 
lines represent the total emission from the disc on a 6 AU field-of-view, respectively with and without the scattered emission. The 
photospheric SED, simulated by a NextGen model atmosphere, is represented by the two solid lines as seen pole-on (upper curve) 
and from the equatorial plane (lower curve). Only measurements plotted in red have been taken into account in the fit. 



up" of the whole planetary system, similar to the late heavy bom- 
bardme nt (LHB) t hat occurred in the early history of the solar 
system ( Gomes et a l. 2005), could enhance this transfer rate and 
induce the necessary amount of dust in the inner planetary sys- 
tem (the solar zodiacal cloud is supp osed to hav e been up to 10"* 
times brighter during the LHB, Nesv ornv et aLl bOlO). Although 
the presence of giant planets around Vega has not been con- 
firmed yet, such a bombardment would most probably be trig- 
gered by_flie out ward migration of g iant planets. In particular, 
IWvattI (120031) and iReche et al] (120081) suggest that the outward 
migration of a Neptune to Saturn-mass body from 40 to 65 AU 
could explain the observed clumpy structure in Vega's outer disc. 
Further dynamical analyses show that such a transport would re- 
quire the presence of at least two planets, such as a Saturn-like 
planet orbiting between 40 and 65 AU and a Jupiter-like planet 
orbiting closer to Vega (Vandeportal et al. 201 1, in preparation). 



6. Conclusion 

Interpreted as the signature of hot exozodiacal dust, the 
near-infrared excess source detected around Vega by 
CHARA/FLUOR in the K band is confirmed by our 
IOTA/IONIC measurements at the 3-cr level. Using the 
new constraints provided by the H-band data (including both 
high-accuracy visibilities and closure phases) and a significantly 
improved spatial coverage, the most straightforward scenario 
consists in a compact dust disc producing a thermal emission 
that is largely dominated by small grains located between 0.1 
and 0.3 AU from Vega and accounting for a relative flux with re- 
spect to the stellar photosphere of 1.23 + 0.45%. This flux ratio 
is shown to vary slightly with the geometry of the model used 
to fit our interferometric data (variations within +0.19%). Using 
the new H-band flux ratio together with archival measurements 
at various wavelengths, we show by means of SED modelling 
that at least a small fraction of carbonaceous grains must be 
present in the disc in order to fit the data in a satisfactory way. A 
steep density profile is also necessary to ensure the compatibility 



with the non-detection reported by nulling interferometry in the 
mid-infrared. Considering the best-fit model, a dust mass of 
approximately 2 x lO^^M®, equivalent to the mass of an asteroid 
about 20 km in diameter, would be necessary to explain the 
near-infrared excess emission. Given the short lifetime of dust 
in the inner Vega system, a major dynamical event, similar to 
the solar system's late heavy bombardment, might be currently 
ongoing in the Vega system. This would support the idea that 
the debris disc around Vega is the only visible component of a 
more complex planetary system harbouring unknown planets. 
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